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A simulation method based on the RG blocking is shown to yield statistical errors 
smaller than that of the crude MC using absolute values of the original measures. 
The new method is particularly suitable to apply to the sign problem of indefinite 
or complex- valued measures. We demonstrate the many advantages of this method 
in the simulation of 2D Ising model with complex- valued temperature. 



1 The Sign Problem 



^ , In order to evaluate a multi-dimensional integral of the partition function 
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O^ . using Monte Carlo (MC) one can sample the points in the integration domain 

^^ ' with a non-uniform distribution, p. This sampling gives the following estimate 

for the integral: 
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pH I where N is the number of points sampled, p > and is normalised; 
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The best choice of p is the one that minimises the standard deviation 
squared S. This can be found by variational method leading to the crude average-sign MC weight 

Pcrudc - \f\/J\f\dV, (4) 
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The sign problemu arises when Z/ J \f\dV is vanishingly small: then un- 
less a huge number of configurations are MC sampled, the large statistical 
fluctuations of the partition function render the measurement meaningless. 



2 Improved Simulation Method 

One way of smoothing out the sign problem is to do part of the integral ana- 
lytically, and the remainder using MC. The analytical summation is not just 
directly over a subset of the dynamical variables; in general it can be a renor- 
malisation group (RG) blocking where coarse-grained variables are introduced. 
This does yield certain improvement over the crude MC in general. 

Let P{V', V} be the normalised RG weight relating the original variables 
V to the blocked variables V, 



P{V',V} > 0, 
I P{V',V}dV' = 1. 



Inserting this unity resolution into the integral (|l|) 

/ = JdvJdV'P{V',V}f, 
dV'giV). 
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Thus, an MC estimator is only needed for the remaining integration over V' 
in m) . As with the crude method of the last section, variational minimisation 
for 5 of (|^), with g in place of /, leads to the improved MC 
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Firstly, the improved weight sampling yields a partition function of mag- 
nitude not less than that sampled by the crude weight: 
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Secondly, it is also not difficult to sec that the statistical fluctuations as- 
sociated with improved MC is not more than that of the crude MC, 
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Thus the RG blocking always reduces the statistical fluctuations of an 
observable measurement by reducing the magnitude of VS/\{{sign))\. 

Note that the special case of equality in (Hpl) occurs iff there was no sign 
problem to begin with. How much improvement one can get out of the new MC 
weight depends on the details of the RG blocking and on the original measure 
/■ 



3 Complex 2D Ising Model 

The Hamiltonian for the Ising model on a square lattice is 
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For the finite 



Here we allow j and h to take on complex values in general, 
lattice, periodic boundary conditions are used. 

The phase boundaries for the complex temperature 2D Ising model with 
/i = are depicted in Figure 1 El. 

For the improved method, we adopt a simple RG blocking over the odd 
sites, labeled o. That is, the analytic summation is done over the configuration 
space spanned by the even o sites; while MC is used to evaluate the sum over 
the remaining lattice of the odd • sites. Figure 2 shows the two sublattices, 
and how the • sites are to be labelled relative to the o sites, for the site labelled 

X. 
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Figure 1: Phase diagram in the complex j plane. 
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Figure 2: Relative spin positions on a partitioned lattice. 
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Figure 3: | << sign >> j vs Re{j) Im(j) = 0, /i = 0. Filled boxes are the improved data. 



In general, with finite-range interactions between the spins, one can always 
subdivide the lattice into sublattices, on each of which the spins arc (nearly) 
independent for the partial sum to be carried out. 

Summing over the spins So, 
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where s^ is the sum of the spins of the four odd neighbours surrounding each 
even site. The improved MC weight is then the absolute value of the summand 
on the right hand side of the last expression for Z . 

The quantities measured are magnetisation, M, and susceptibility, x- 
These can be expressed in terms of the first and second derivatives of Z re- 
spectively, evaluated a,t h — Q. In all the simulations, square two-dimensional 
lattices of various sizes with periodic boundary conditions are used. The 
heat-bath algorithm is used. Two additional benefits arise from the improved 
method. The first is that the number of sites to be visited is halved. While 
the expressions to be calculated at the remaining sites turn out to be far more 
complicated, the use of table look-up means that evaluating them need not 
be computationally more expensive. The second benefit is that the number of 
sweeps required to decorrelate data points is reduced. 

As a test of the improved method, it is compared to the crude one along 
the path OX in Figure 1. Comparison with series-expansion data are plotted 
in Figure 4, in which the discs are the absolute- value circles of complex- valued 
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Figure 4: Improved x vs Re{j): expansion data shown as line. 

statistical errors for simulated results. 

4 Concluding Remarks 

The sign problem is partially eleviated in the improved method H because of 
some partial phase cancellation among the original indefinite or complex- valued 
measure after an exact RG transformation. A particular RG blocking is chosen 
for our illustrative example of the 2D Ising model with complex- valued mea- 
sure. But other choices of RG blocking are feasible and how effective they are 
depends on the physics of the problems. 
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